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We show that zero ejection energy of the photoelectrons is classically impossible for hydrogen-like 
ions, even when field ionization occurs adiabatically. To prove this we transform the basic equations 
to those describing two 2D anharmonic oscillators. The same method yields an alternative way to 
derive the anomalous critical field of hydrogen-like ions. The analytical results are confirmed and 
illustrated by numerical simulations. 
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I. INTRODUCTION 



O 

A wealth of new phenomena in the interaction of strong, short laser pulses with matter has been found in recent 
years Laser systems delivering pulses at irradiances up to I0 19 Wcm~ 2 ^m 2 are available in several laboratories. 

Pulses not longer than a few tens of femtoseconds are "state of the art" . In describing their interaction with matter 
usual perturbation theory fails since the electric field of such laser pulses is of the order of the atomic field. One 
^ • of the prominent phenomena occuring in a strong laser pulse is field ionization. Here the electron is able to escape 
over the barrier formed by the Coulomb potential and the laser field. This process is much faster than tunneling 
which is dominant in weaker fields. In other words, in intense laser fields the so-called barrier suppression ionization 
(BSI) regime jj) is reached. In this regime, classical and semi-classical pictures are expected to work well and much 
. theoretical research has been done in that direction, both analytically Q-j^] and numerically ]7|-|ri]]. 
' So-called "simple man 's theory" (SMT) 0-0 succeeded in explaining essential features of above threshold ion- 
ON '. ization (ATI) (see e.g. |j],[n| for an overview), such as cut-offs in photoelectron energies and harmonic spectra Jl6| . 

■ Extensions to SMT considering rescattering effects clarified experimentally observed features such as plateaus in 
O \ ATI-spectra |l7]] and rings in the angular distributions p8j . 

* 55 ■ However, no ionization dynamics enters in the SMT since the electron is regarded as being "born" in the laser 
field without perceiving any attraction by the ionic core anymore, at least until rescattering by the nucleus. Even in 
, Keldysh-Faisal-Reiss (KFR)-type theories Jl^ the Coulomb potential enters not directly but through the initial 
Oh- energy state wave function. 

For relevant applications, e.g. inverse bremsstrahlung, it is important to know the energy with which the electron 
is freed (ejection energy) p2j. To calculate the latter the ionization dynamics must be described by solving the 
time-dependent Schrodinger equation or by using simpler alternative concepts. In this paper a classical model is used 

■ for this purpose. 

For the sake of completeness, we briefly review SMT. Suppose that a linearly polarized laser field E(t) sets in at 
t = 0, the velocity of a non-relativistic electron released at time t' is at the instant t 

«(*) = -/ E(t) dr + v o = - I E(T)dT + v {t')+v Q (I) 



o 



(atomic units will be used throughout this paper) with 

,o 

«d(0 = -J E(r)dT 
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the residual drift due to ejection out of phase, i.e. not in the maxima of the field. i>oo is the ejection velocity we are 
particularly interested in. Since the ejection velocity vqq contributes to the overall residual kinetic energy £ rcs of the 
electron at the end of the pulse, 

£rcs = ^m(v + Woo) 2 , 

it alters the electron distribution function and must be taken into account. Ignoring vqq in (|l|) leads in the case 
of linear polarized laser light to the well-known results (£ re s) = 3C/ p , max£ re s = 8U P as well as to the maximum 
return energy w 3.2U p which is responsible for the cutoff in harmonic-spectra. U p is the ponderomotive potential 
e 2 E 2 / (Amu> 2 ) , i.e. the mean quiver energy of the electron in the field. 

The paper is organized as follows: in Sec. [n| the anharmonic oscillator picture is introduced which leads to a 



physically appealing interpretation of the adiabatic ionization process as discussed in Sec. [II . In Sec. IV, the lower 
limit for the ejection energy as well as the upper limit for the ejection radius is derived. It is further shown that no 
upper limit for the ejection energy exists. In Sec. [v|the analytical results are confirmed and illustrated by numerical 
simulations. Sec. [v| is devoted to a rederivation of the anomalous critical field of hydrogen like ions |2(J in the 
framework of the anharmonic oscillator picture. The adiabatic deformation of the initial ground state is calculated 
and compared with the numerical solution of the time-dependent Schrodinger equation. Finally, we give our conclusion 



in Sec. VII 



II. TRANSFORMING TO ANHARMONIC OSCILLATORS 

The laser frequency Q can be regarded as small compared to the orbital frequency uj of the bound electron (e.g. 
il = 0.18 a.u. for KrF and £1 = 0.04 a.u. for Nd). Therefore, a quasi adiabatic treatment is appropriate. We will now 
outline the transformation to a set of two 2D anharmonic oscillators 



We start with the Hamiltonian describing an electron in a Coulomb potential —Z/y/p 2 + z 2 , and a static electric 
field E aligned in z direction, 

H(p,z;p p ,p z , P(p ) = ± (pl+pl + ^j -—A— + Ez = E. (2) 

Here, £ is the total energy of this conservative system and the azimuthal momentum p v is a constant of the motion. 
We transform to parabolic coordinates first. The new coordinates £ and r] are related to the cylindric ones according 

to 

£=(r-z)/2, V =(r + z)/2, 
with r = \J ' p 2 + z 2 . The canonical momenta are related through 

This leads to the Hamiltonian in parabolic coordinates, 

1 ( i 2 , V 2 , 1 2\ z 




It is well known that the Hamilton- Jacobi-Equation of the problem separates in parabolic coordinates Jm) , as does so 
Schrodinger's equation. It is advantageous to perform another canonical transformation p3| , 

£ = u 2 /4, 77 = w 2 /4, < u,v < 00, 

Pu = \/lpe,, Pv = y/rjPr,, 



leading to 

H{u,v;p u ,p v , Pip ) = - ( 2 ^ 2 J [pI + ^pI +Pl + ^pI ) - 7^-71 + 9{v 2 - u 2 ) = £, (4) 
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where the new field g = E/4 has been introduced. 

Now, we define the "zero energy Hamiltonian" Ho = H—£ = 0, and by multiplying this equation with (u 2 +-u 2 )/4 = r 
we are left with 

H o = \(pl + §+p" + P ^)+Y {u2 + v2) - z + \ 9{v4 ~ ui) ^^ 



where uj 2 = —£/2 . Therefore, we finally get 

if p 2 \ uj 2 1 

H u (u;p u , Ptp ) = - I pi + -§ J + —u 2 - -gu A = A, (5) 

H v (v;p v , P(p ,) = ~ + + ^_ v 2 + ]_ gv i = B (6) 

A + B = Z, (7) 

P V =P V >, (8) 



which represents a set of two 2D anharmonic oscillators moving independently. There are only constraints concerning 
the initial values of the system: (i) the energies must sum up to the given total "energy" and (ii) the angular 
momenta, which are constants of the motion, are equal. 

III. PHYSICAL PICTURE 

The physical interpretation of Eqs.(|^) and (^J) is very simple. Supposing g > 0, Eq.(^|) allows only for a bound 
motion in v, whereas the effective potential in Eq.(||), V u (u) = p 2 /(2u 2 ) + to 2 u 2 /2 — gifi/A, has a local maximum so 
that the motion in u is unbound if the energy A lies above the potential barrier. 

Now, the field g may be turned on adiabatically. Initially, when g — holds, the v- and the u- motion take place in 
equal potentials. While g increases, the u-potential is bent down whereas the w-potential is steepened. Therefore one 
expects an adiabatic lowering of the energy A and an adiabatic raising of the level B (Fig. |l|). But the frequency lo 
may also change because the adiabatic invariant 

S u {uj,g,A) = V2 <j> V 'A - V u {u)du = const. (9) 

yields a relation between g, A and to only. Indeed, a change in to is the classical dc Stark-shift, of course. The trick 
is to consider the special Kepler orbit which ionizes earliest in an adiabatically ramped dc field. Therefore, we switch 
back to the "physical space" for a moment: the classical trajectory which ionizes first is that one mostly directed 
towards the potential barrier. There are two Kepler orbits lying on the z-axis in the limit eccentricity — > 1. One lies 
in the region z > 0, the other one in z < 0. If the potential barrier comes from negative z, as in our case where g > 0, 
the Kepler orbit z < ionizes first and the Kepler orbit z > does last. The two z-directed orbits are the extreme 
Stark-shifted ones. The orbit which ionizes first corresponds to the electron with smallest ejection energy. Therefore, 
the ejection energy of just this electron provides a lower limit for all ejection energies which may occur. Fortunately, 
this special energy can be calculated. 

Now we turn back to the anharmonic oscillator picture: motion along the negative z-axis means v — p v — p v > = 
for all times. Therefore, the anharmonic oscillator H v is frozen, i.e., B = 0. Now, when g is adiabatically ramped B 
is not adiabatically raised since the v-particle rests in the potential well and does not recognize the steepening of the 
potential. It follows that B remains zero, and A = Z (Fig. ||). The adiabatic invariant (||) now provides a relation 
between u) and g only, i.e., one can, in principal, calculate the classical Stark-shift 0J 2 {g) for this particular orbit. The 
Stark-shift is the key which enables us to calculate the lower limit for the ejection energy and the upper limit for the 
ejection radius, as we will see in the following Section. 
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IV. CALCULATION OF THE EJECTION ENERGY 



Since the integral in Eq.Q can not be solved analytically some approximations must be applied. We will restrict 
ourselves to the linear Stark-shift which follows from Eq. (S) when the integrand and the turning point is expanded in 
powers of the field g. The result corresponds to the quantum mechanical one in the limit of high quantum numbers, 
n — ► oo. 

The potential barrier for H u in the case p v = is located at u& = u>/ \fg. The value of the potential at that point 
is V u (ub) — ut 4 / (Ag). Since the kinetic energy is zero when the potential barrier meets the total energy A = Z for the 
first time (this is the case due to adiabaticity) , one has 

%=A = Z (10) 

at that certain moment. 

Now, we assume a linear dependence between energy and field, i.e., only the linear Stark effect is included in our 
further calculations, 

£(E)=£ + ^E, A<0. (11) 
A 

As derived by expanding the invariant (^) or from the classical limit of the well-known quantum mechanical result 

11 
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holds, which will be used in the following. 
The "over the top" -criterion (I1Q) then reads 



and this leads to the "over the top" -field strength 

E b = E(£ b ) = ^. (12) 

The electron is not yet free when it flows over the barrier because its Coulomb energy £ — Ez is still negative. 
Demanding £ — Ez = 0, which may be written in u and v as £ — g(v 2 — u 2 ) — 0, is equivalent to claim 

V u {u) = ^u 2 - ^gu 4 = 

when v = 0, as it is in our case. Therefore = A = Z holds during the moment when the electron becomes free. 
The coordinate Uo where this happens is easily determined through Uq/A = uj 2 /(2g) This equation gives an upper 
limit for the ejection radius, 

ro = 3^. (13) 

The physical kinetic energy T is related to p u according Eq.(|4|) (using p v = v = p v = 0) through T := 2p 2 l /u 2 . 
For uo this leads to To = 2Zg/u 2 . Going back to the physical quantities £ and E instead of to 2 and g this means 
To = —ZE/£. Assuming that £ and E do not vary much during the short time between "flowing over the top" and 
"getting positive Coulomb energy" , i.e. ionization, and inserting A, we get 

To = \\£o\- (14) 

Note that the ejection energy is essentially the kinetic energy gained by the "u-particle" when it slides from the top 
of the potential barrier down to V u (uo) = 0. Multiplication with a factor 4/mq leads to the physical ejection energy 
T . 
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We have already mentioned the two extreme type of orbits: they are both aligned along the z-axis (eccentricity 
— > 1), the early ionizing one on the negative half, the latest ionizing orbit on the positive half. For the latter one 
u = and A — holds instead of v — and B — 0. Therefore, we have only to substitute g — > — g in the result for 
the Stark-shift. Obviously, there is no potential barrier in V v (v) over which the electron could escape. Nevertheless, 
from the "zero Coulomb energy" -condition V v (v) = follows, as in the it-case. This condition can only be fulfilled 
through v = 0. Thus, u = v = holds simultaneously, i.e., the electron ionizes over the Coulomb singularity in the 
limit eccentricity — > 1, and the kinetic energy is infinite at that point. So we conclude that there is no upper limit for 
the ejection energy. 



V. NUMERICAL RESULTS 



Our numerical simulations presented in the following were performed for the "classical" Is state of atomic hydrogen, 
Z = 1, £q = — 1/2. The classical linear Stark effect then gives A = —8/3. The "over the top" -values are calculated to 
be £b = —2/3, Ef, = 1/9, and, finally, the lower limit for the ejection energy is found to be To = 1/6. Eq.([l3|) gives 
for the ejection radius ro = \zq\ = 6, i.e., all electrons will be ejected within a sphere of radius 6 a.u.. 

Numerically, the exact motion of an ensemble of electrons in a Coulomb potential and an electric field can easily be 
determined using the so-called "classical trajectory Monte Carlo-method" (CTMC) [^5|. We have performed CTMC 
runs and looked at each test electron when its Coulomb energy £ — E(t)z becomes positive, i.e., when the electron 
would be free if the field is turned off immediately (this is our definition of ionization) . 

In Fig.|| the kinetic ejection energy of each electron is indicated by a symbol over the field strength where ionization 
takes place (the field is directly proportional to the time due to linear ramping) for Z = 1, £q = —1/2, E(t) = t/2000, 
and an ensemble consisting of 1000 test electrons. The lower limit of the ejection energy is indicated in the plot by a 
horizontal line. 

In Fig]!] a similar plot for the ejection radius is presented. 

An oscillating field E(t) ~ cosftt can be treated adiabatically when the frequency £1 is small compared to the 
orbital frequency |£o|/?i, i-e. in atomic units O / 1 f 1 <C 1. We repeated our numerical run discussed in the previous 
paragraph with the field replaced by E(t) — t cos(fit)/2000, ft — 0.04. The frequency ft corresponds to Nd laser-light 
and due to f2/|£o| = 0.08 an almost adiabatic behavior is expected. In Figj| the kinetic ejection energies of each test 
particle is plotted vs. its ionization time. The absolute value of the electric field and the calculated lower limit for the 
ejection energy are included. Apart from three test particles which are released extremely out of phase, i.e., when the 
electric field is low, all electrons lie above the calculated lower limit. 

VI. QUANTUM MECHANICAL CALCULATIONS 

The anharmonic oscillator method provides also an elegant and alternative way to derive the anomalous critical field 
p6| for hydrogen-like ions. With the ansatz ip(u, v, tp, ip'(ip)) — ipu(u) exp{irrup)ijj v {v) exp(imV') the Hamiltonians (|^) 
and (H) lead to the following two Schrodinger equations 



-d u (ud u ) - ^ 



-u - -gu > ip u (u) = Aij} u (u), 



-d v (vd v ) - 

V V 
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m = m , 



A + B = Z. 



Since the solution of the unperturbed problem, which is a 2D harmonic oscillator, is a Gaussian in the m = m! = 0- 
case, we use 



ipu(u) = V au/TT exp(-a u u 2 /2) 



as a trial function. The total energy then is 
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(I 



M.!l) =2,t / .!,«.,;/[„, „ = J-^-M +««-^2 



Minimizing this energy yields up to first order in g a u = uj(1 — g/u 3 ) and a v = u>(l + g/uj 3 ). The oscillator energies 
are A(g) = lo(1 — g/(2u> 3 )) and B(g) = u(l + g/(2co 3 )). Note that this is consistent with the fact that the linear 
Stark-effect vanishes for the ground state of hydrogen-like ions, since A(g) + B(g) = 2lu and from A + B = Z and 
to 2 = —£o/2 follows £q — — Z 2 /2. For the hydrogen ls-state exposed to an adiabatically ramped dc field E, the wave 
function in physical coordinates is 



1-4J5 
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i/> H (p, z) = ^— exp(- Vp 2 + z 2 ) exp(-2£z), (15) 

i.e., the unperturbed wave function is multiplied by a "deformation factor" exp(— 2Ez). 

In order to calculate the critical field £? C rit = 4<? cr it we claim that the u-Gaussian touches the potential barrier 
energetically: 

4 

A(g) = lu- -2L = . 
y ' 2uj 2 Ag 

This leads to 

5crit = L0 3 (1 - 1/V2) P3 0.3w 3 . 



For the Is hydrogen state we get E^ it = 0.15 a.u. in accordance with the result presented in p6|] . 

In Fig. ^ a comparison is made between the calculated deformation \iPh(p, z )\ 2 which results from Eq. [l5] and 
the exact numerical solution of the time-dependent Schrodinger equation where we ramped the field linearly over 30 
atomic time units up to the critical field E^. it = 0.15 a.u.. Strong deviation occurs only for the 1%-contour while the 
main parts of the probability densities agree well. 

In Fig. the initial probability density and the variationally determined one at the critical field are shown in the 
energy diagram. Note that at the critical field the packet lies 0.27 atomic energy units above the barrier in z-direction. 
This value may be interpreted as a quantum mechanical ejection energy. 



VII. CONCLUSIONS 



We have shown that within the framework of the anharmonic oscillator model important features in strong field 
ionization can be derived. Using this method, we have calculated (i) the lower limit for the ejection energy, To = |£q|/3, 
(ii) the upper limit for the ejection radius, ro = 3Z/\£q\, and (iii) the anomalous critical field for hydrogen-like ions 
as well as the shape of the deformed probability density. We have also demonstrated that there exists no upper limit 
for the ejection energy. All results have been confirmed by CTMC simulations. 

FIG. 1. 

The two 2D anharmonic oscillators in the case p v = p v > = 0. Initially, when g = holds the two potentials are 
pure harmonic ones (dashed). As g is adiabatically raised one expects an adiabatic lowering of level A and a raising 
of level B. 



FIG. 2. 

For the early ionizing Kepler orbit discussed in the text v = holds. Therefore, and owing to the constraint 
A + B = Z, both levels do not move when g is adiabatically raised. 



FIG. 3. 

The ejection energy of all test electrons in z- and p direction (O, A), and the sum of both (■). The calculated lower 
limit for the ejection energy is also shown. One can see that this limit is confirmed by the numerical runs: none of 
the kinetic energies ■ lies beneath the limit. 
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FIG. 4. 

The radii of all test electron when they become free. The analytically calculated upper limit is confirmed by the 
numerical run. All electrons are ionized within a sphere of radius 6 atomic units. 



FIG. 5. 

The kinetic ejection energies in an oscillatory i?-ficld. Only 3 of 1000 test electrons ionize through a "non-adiabatic 
channel" in such a way that they lie beneath the calculated lower limit. Note, that these particles become free when 
the field is low. The absolute value of the electric field is also shown. 



FIG. 6. 

The contour lines representing 80, 50, 10 and 1% of the peak height at the origin p = z = are shown for the 
numerical result (solid) and the variationally calculated probability density (dotted). The field strength is the critical 
one E = 0.15 a.u.. 



FIG. 7. 

Unperturbed and Stark-deformed Is peak for atomic hydrogen at the critical field strength E CTl t = 0.15 a.u.. The 
energy gap of 0.27 a.u. between the energy level at £q = —1/2 and the barrier can be interpreted as an quantum 
mechanical ejection energy. 
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